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Abstract-—Pour-compartment models in the form of a system of linear differential equations, 
associated with in vivo cadmium data, have been constructed to describe a particular biological 
phenomenon of cadmium metabolism in the human body. Since it is quite technically expensive 
to make in vivo measurements and also difficult to collect data over a period of a decade, the data 
analyzed in this paper are quite precious, and in fact, one of very few such existing data sets in the area 
of cadmium research. Cadmium researchers are very much interested in drawing as much information 
as possible from these data. Methods are developed here for deriving the expectation functions 
and for using them to analyze this special data set. The parameter estimation for compartment 
models discussed next uses classical as well as Bayesian approaches. This is the first time that a 
whole system of the human body is analyzed simultaneously, instead of discussing each compartment 
separately with a large number of assumptions from different sources on Kjellstrom's model, without 
any additional assumption on the derived compartment models. The results, obtained from these 
statistical approaches on simple as well as more direct mathematical models, are not only interpretable 
very reasonably in terms of biological phenomena, but also reveal interesting patterns and show a 
great consistency with previous studies. © 2000 Elsevier Science Ltd. All rights reserved. 


Keywords —Cadmium metabolism, Compartment model, Bayesian estimation, Parameter esti¬ 
mation. 


1. INTRODUCTION 

Cadmium has been used in a wide range of industrial applications even though its toxic ef¬ 
fects have been acknowledged for decades [1]. Toxic effects of cadmium occur if some significant 
biological mechanism is stimulated beyond a certain level, Of special interest among chronic 
effects of cadmium is kidney damage (being the most typical feature of chronic cadmium intox¬ 
ication), which is characterized by an excess excretion of low molecular weight proteins such as 
/? 2 -microglobulin [2]. 

Cadmium researchers are therefore quite interested in biological phenomenon of cadmium 
metabolism, and in particular with the absorption, distribution, accumulation, and excretion 
of cadmium with respect to time in the human body. The absorption and subsequent distribu¬ 
tion of cadmium depend on the particular exposure route and the form and intensity to which 
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workers are exposed. Cadmium is usually inhaled or ingested and, initially, a considerable amount 
will be found in lungs or intestinal wall from where it is transported, by blood plasma, to other 
tissues. Cadmi um in bloo d plasma is bound to low molecular weight as well as high molecular 
weight proteins, and it is the transport and metabolism of these that determines, to a large 
extent, the distribution of cadmium through the body. The main cadmium accumulation or¬ 
gans are the liver and kidneys, which account for approximately 16% and 53% of the total body 
burden, respectively [2,3]. Accumulation in these sites results from the binding of cadmium to 
proteins—particularly, albumin and metaliothionein for liver, and raetallothionein for the kidney. 
An eight-compartment kinetic model, showing a flow scheme of cadmium metabolism, has been 
established. This model, with 21 distribution coefficients, contains compartments of liver, kid¬ 
ney, feces, urine, other tissues, b lood 1, b lood 2, and blood 3 cadmium. Blood 1 contributes to 
accumulation of cadmium in liver and other tissues, while blood 3 contributes to accumulation 
of cadmium in kidneys. Blood 2, between blood 1 and blood 3, contains the accumulation of 
cadmium bound to cells and molecules. Background information for Kjellstrdm’s model had been 
obtained partly from animal experiments and partly from observations on the human beings in 
industrial and general environmental exposure situations. 

In addition to a number of assumptions in Kjellstrom’s model, all the unknown parameters are 
fixed according to prior experience and knowledge based on different sources and different animal 
and autopsy data sets. During the analysis, the numerical solution of only one compartment of 
the model can be derived after the unknown parameters are fixed. Thus, the relationship between 
each compartment and time can be established with prior fixed parameters. In contrast to Kjell¬ 
strdm’s model, the four-compartment model discussed in this paper is formulated with systems 
of linear differential equations in a simple and more direct manner. The parameter estimation 
for the four-compartment model is discussed using classical and Bayesian approaches. This is 
the first time that a whole system of the human body is analyzed simultaneously without any 
additional assumption on the derived compartment models; usually, discussions are made on each 
compartment separately, with a number of assumptions concerning the unknown parameters. Fur¬ 
thermore, all the parameters are estimated from observed human data in the four-compartment 
model. 

Kjelistrom’s model provides a good biological basis and reference for the starting values of the 
unknown parameters while dealing with the four-compartment model, even though the parameter 
estimation based on the four-compartment model is achieved in an inverse way mathematically. 
The results obtained here are not only interpretable very reasonably in terms of biological phe¬ 
nomena, but also reveal interesting patterns and show a great consistency with previous studies, 
although the fact that data here are sparse and noisy brings about the special difficulty encoun¬ 
tered in this study. 


2. PROBLEM FORMULATION 

2.1. Subjects 

The subjects of this paper are 14 male workers involved in brazing using cadmium-silver alloys. 
The duration of cadmium exposure ranged from 2 to 34 years. The level of exposure for each 
subject varied throughout the exposure time according to the work undertaken. The exposure for 
all but one subject (number 3) was at or above the current U.K. occupational limit for cadmium 
in air, i.e., O.G5mg/m a , for at least part of the time. 

Longitudinal measurements including biochemical and in vivo measurements were made for 
these 14 male workers. Biochemical measurements consisting of urinary cadmium, b lood cad ¬ 
mium. and urinary f^-microglobulin had been made at approximately six monthly to yearly 
intervals from 1983 to 1990. Blood and urinary cadmium were measured by using electrothermal 
.atomic spectrophotometr y, while the urinary /^-microglobulin was measured in freshly voided" 
samples using a radioimmunoassay procedure. 
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In addition to the biochemical measurements, all subjects had in vivo measurements of liver 
and kidney cadmium made in December 1983, and eight subjects (numbers 1 , 2 , 3, 4, 5, 10, 
11, 13) had additional in vivo measurements made in March/April 1990. These organ levels were 
measured by prompt 7 -ray neutron activation analysis, using different measurement systems on 
each occasion. The differences arising solely from the different measurement systems are thought 
to be negligible. 

In vivo measurement of the levels of liver and kidney cadmium along with urine and blood 
cadmium, as well as low molecular weight proteins /? 2 -microglobulin provide estimates in order to 
describe the phenomenon of cadmium metabolism, based on the compartment mode! [4). Since 
it is quite expensive to make in vivo measurements and also difficult to collect data over a period 
of decades, the data analyzed in this paper are quite precious, and in fact one of very few such 
existing data sets in the area of cadmium research. Cadmium researchers are interested in drawing 
as much information as possible from thesejlgta. 

In this paper, the abbreviations Cd_u, C’d-b) Cd_l, CdJr, and uJ ?2 are used in some formulae, 
which denote measurements of urine cadmium, blood cadmium, liver cadmium, kidney cadmium, 
and urinary /Jg-microglobului, respectively. The duration of cadmium exposure, the first mea¬ 
surement and last measurement of u ./?2 as an indicator of kidney status, the hi story of smokin g 
and age for each subject are listed in Thble 1. 

Table 1. Duration of exposure, age, and indicator of kidney status. 


Subject 

Number 

Cadmium 

Exposure 

(years) 

Year of 
Birth 

*/? 2 -MicroglobuIin 
(mg/g ercat.) 

Smoking Habit 

1 

2 

1964 

257/NA 

not smoking 

2 

31 

1937 

90/57 

not smoking 

3 

13 

1945 

58/135 

not smoking 

4 

33 

1930 

86/238 

smoking (1948 -) 

5 

20 . 

1928 

688/88 

not smoking 

6 

9 

1923 

344/NA 

NA 

7 

3 

1961 

19/NA 

NA 

8 

2 

1925 

59/NA 

NA 

9 

24 

1928 

1667/NA 

NA 

10 

34 

1921 

3077/84600 

smoking (1939-1958) 

11 

16 

1920 

158000/154100 

smoking (1938-1962) 

12 

22 

1927 

13400/155800 

NA 

13 

22 

1925 

1067/NA 

smoking (1941-1975) 

14 

24 

1924 

37000/102500 

NA 


In the ‘/32-microglobulin’ column, the left side of backslash is for the first 
measurement, and the right side is for the last one. 

‘NA’ means that the information is not available. 


2.2. Compartment Model 

Derivation of the solution of the four-compartment model, a linear system of ordinary differen¬ 
tial equations, is discussed in this section, in which the four compartments are ca dmium in blo od, 
cadmium in urine, cadmium in kidney, and cadmium in liver. It is assumed that the rates of 
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flow of cadmium follow the first-order kinetics, in which the mass balance equations are required 
for the description of dynamics of the exchange of cadmium among compartments. The transfer 
coefficient is assumed to be constant with respect to time. The model can then be described 
analytically in the form of a linear system of ordinary differential equations as follows: 


= AH + B, t> 0, 
or 

H( 0 ) = tfo. 


where 


H = 


'%■ 


— (Cm + C&*) 

Cib 

C hb 

O' 

Vi 

A = 

Cu 

"Cib 

0 

0 

Vk 

i — 

Cbk 

0 

— (Cfcb + Cfcu) 

0 

.n*. 


0 

0 

Cku 

0, 


-R' 

0 

0 ’ 
0 


with the parameters Cm, Gy,, Cbk, Ckb> and C* u in the matrix A denoting the transfer rates 
of cadmium from blood to liver, from liver to blood, from blood to kidney, from kidney to 
blood, and from kidney to urine, respectively; R in the vector B stands for the monthly intake 
of cadmium absorbed into blood from environment; the elements of H, viz. tjb, ift, 17 %, and tj u , 
denote the expected responses of cadmium in blood, liver, kidney, and urine, respectively; t is 
time in monthly unit; and each element of Ho denotes the initial value of corresponding element 
of H at time 0. 

The compartment system can also be represented graphically, in which the boxes represent 
compartments, and the arrows labeled with transfer coefficients represent transfers of material 
into or out of compartments. The compartment diagram for the model described above is shown 
in Figure 1. 


i Cd 
| Blood 

4-t 

Liver 

_ ^ 
Urine 


n 

Kidneys 


Figure 1. Diagram for the four-compartment model (Cd.b, CdJ, Cd-k, and CtLu). 


Since the measurements on the amount of cadmium in urine per month are available and urine 
cadmium has little feedback to kidney, the dimension can be reduced to take the first three 
ordinary differential equations into consideration at the initial 6tage of deriving equations for the 
expected responses. For simplicity in computation, we shall use 


where 


dHi 


— Ai Hi + Bi, 


dt 
dvu 


dt 


= CkuVk, 


t> 0, 


'Vb' 


'~(Cu + C b k) 

Ci g 

C kb 


~R' 

m 

, j4i — 

C bl 

~C,b 

0 

, j 

0 

.Vk. 


. c bk 

0 

— (Ckb + Cku). 
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The common approach for constructing the expectation functions is to derive the weighted 
sum-of-exponential model in which the expected responses can be expressed as combinations of 
■weighted sum of exponentials and the constants due to the special form of the ordinary differ¬ 
ential equations (without solving them). From this approach, one cannot directly estimate the 
parameters in the generating system, and the number of exponential parameters may exceed the 
number of system parameters. It is not practical to use this approach for setting the functions 
of expected responses on the cadmium data we have here, since the number of measurements for 
each subject is at most 12 , which is not large enough, and the data structure is irregular with 
some missing values. 

We, therefore, adopt here a combination of analytical approach and exponential matrix ap¬ 
proach. Under this approach, the expectation functions for the blood, liver, and kidney cadmium 
can be expressed by elements of A\, the elements of Bi, the eigenvalues of Ai, and the elements 
of the corresponding eigenvectors. 

In order to find all the eigenvalues of A\ and the corresponding eigenvectors, consider the 
characteristic polynomial of Ai given by 

p(A) = deb (Ai - A/ 3 ) , 

where det(j4i — A/ 3 ) denotes the determinant of matrix Ai — A/ 3 , I 3 denotes a 3 X 3 identity 
matrix, and A is a scalar. This characteristic polynomial can be simplified as 

A 3 4- mi A 2 + m 2 A + m 3 — 0 , 


where 


mi ~ Cbt 4- Cu, 4 - Cbk + C k b + Ckv., 

m 2 = C u C kb + CbiCku + C bk C kb 4- C'tfcCjiu 4- CibCku + Ci&C&fc, 
m 3 = C lb C bk C ku . 


To solve the equation easily, we make the change A = x — mi/3. It then turns out to be 

x 3 4- fix 4- h ~ 0, 

where 

, m? 
f\~ —~ 4- m z , 

2 2 mnn2 , 

/2 = ^m 1 -- r -+m 3 . 


By applying Cardan’s formula [5], the discriminant A for the equation is given by 
There are four types of eigenvalues of Ai, depending on what the above discriminant A is. 


Type 1 

If A is greater than zero, only one eigenvalue is real and the other two are complex conjugate 
values. The expressions of eigenvalues in this case are given by 

A 1 = a4^-y, 

A2 = (a + P) + ^i(a -£) - 
A 3 = (a4j3)-y i{a - /?) - 

and _ _ 

c-^-A + vS; p-yiA-VS, 

where i — \f —T is the imaginary number. 
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Type 2 

If A is equal to zero, one of the possibilities is that /i = h = 0. All the eigenvalues are then 
the same and real in this case, and are given by 

A t = A 2 = A 3 = 


Type 3 

If A is equal to zero, another possibility is that /{* /9 = —/f/4. In this case, two eigenvalues 
are the same and all of them are real. They are given by 


Ai 



and 


A 2 — A 3 


ja + P) _ mi 
2 3 ‘ 


Type 4 

If A is less than zero, the eigenvalues are all distinct. These eigenvalues and the corresponding 
eigenvectors are given by 


and 


where 


-■/Hs-WD 


4- 2 jir 


U 


1 

nj 


1 

n 2 


3 = 1,2,3, 


1 

«3 


\A -(- nf +'pf y/1 + nf + pi y/l+nf+$ 
Pi p2 P3 

1/1 + nf+Pi 1/1 + n| + P 2 \A+ n !+P3-f 

Obi _ , _ Cjfc 


A< 4- Cn, 


and p, = 


Ai 4- Ckb 4- Ok u ' 


i = 1,2,3. 


Biologically, it will not be meaningful to have Type I since it will not be meaningful for 
the eigenvalues, associated with the half-life time of cadmium in blood, kidney, and liver, to 
be complex values. We can rule out Type 2 as well, since it is not realistic to have all three 
eigenvalues be the same, as the half-life time of blood cadmium cannot be the same as the half- 
life time of either liver cadmium or kidney cadmium. There is a diminutive possibility of having 
Type 3 due to the uncertainty in the parameter estimation and the possible similarity between 
the half-life time of kidney cadmium and that of liver cadmium. Therefore, it is most likely to 
be of Type 4 for all subjects in that three eigenvalues derived from matrix A x are all distinct. 

The solutions of linear systems of ordinary differential equations with initial values exist in 
the defined region and can be formed in terms of matrix exponential defined as exp(Ait) = 
see [6]. The matrix exponential can be written as 


exp (Ait) = If exp(Jt)U 1 , 


where Ai (i = are the eigenvalues of Ai, U is nonsingular, the columns of U are the 

corresponding eigenvectors of Ai, and exp(Jt) is diag{exp(Ajt),..., exp(A„t)}, 
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It; is assumed that the three eigenvalues A* (i = 1,2,3) of A 1 are dissimilar in each subject. 
The solutions of the compartment model as expectation functions can then be derived as 


H\(t) = exp (A\t) H \ 0 + exp (At) / exp 1 [An]} B\ dr} 

Jo 

= U exp(J£)17 -1 ifi a 4- Uexp{Jt)U~ i (I — Uexp(~Jt)U~ 1 ) A^Bi 
= UQxp(Jt)U~ l (jfflo + A^Bi) - A^Bu 


t> 0, 

t <0. 

The columns of V are the corresponding eigenvectors of A%. The eigenvalues of A,, A f 
(i = 1,2,3), and the corresponding eigenvectors in the expectation functions which have been 
presented all correspond to Type 4. 

If the situation in Type 3 takes place, the matrix diag{exp(Ait),..., exp(A n t)} should be re¬ 
placed by a general form, viz. the Jordan canonical form. The corresponding eigenvectors should 
be the generalized eigenvectors. This situation, however, will not be discussed in this paper. 

Comparing the expectation functions to the ordinary differential equations, it can be noticed 
that the expectation functions have been modified when the time is less than 0. There are two 
possible reasons for the modification of the model. 

1. Since the cadmium exposure to the workers ceased at time t 0 and the first measurement 
was made when the cadmium exposure already started, the level of cadmium in blood 
cannot be a dramatic change from the time of first measurement to time 0, and it would 
decrease right after the exposure stopped until it reaches another steady state. 

2. There is no direct relationship between the intake of cadmium and liver cadmium, or 
kidney cadmium. The intake of cadmium has direct contact with blood cadmium through 
respiratory system and blood circulation system. 


where 


exp(J"f) — 


0 
0 


0 

e Aat 0 


0 




and 


=r 

lo, 


3. PARAMETER ESTIMATION 


3.1. Two Methods of Estimation 

Multivariate analysis has been used to estimate the unknown parameters with four responses 
(cadmium in blood, liver, kidney, and urine). 

In general, the data for multiresponse parameter estimation consists of M observed responses, 
y nl , y n 2 ,..., y n Mi taken at each of N experimental conditions, x n (n = 1,2,..., N). These data 
are modeled as 

Vnm = Vm (®rti (?) T ^nmi ri = 1, ■ •»i A) m “ 1,..., M , 

where 8 is a P x 1 vector of unknown parameters; all observations y nm are gathered into an N x M 
matrix Y\ ?? m is the model function of known form for the m til response depending on some or 
all of the experimental setting x n and on some or all of the parameters 8 , and all ri m (x n , 8) are 
collected into JV X M expected response matrix H(8)\ e nm is the disturbance term, which is 
assumed to be normally distributed and all the residuals are collected into an N x M matrix 
Z(8) = Y{8) - H{6). 

In the cadmium data, M is always 4, the number of unknown parameters P is 6, N varies 
from 9 to 12, and 8 is the vector [Cu,Cib,Gbk,Ckb,Gk U iW‘ 

Procedures for the determination of parameters in this multivariate system are based on two 
criteria, one based on the generalized nonlinear least-squares minimization, and the other based 
on the Bayesian determinant minimization. 
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The objective function of generalized nonlinear least-squares can be written as 

M N 

h{ 0 ) = "y 1 y ] \ifnm ~ Vm (^rn 0 )] • ( 3 - 1 ) 

m—1 n=l 

The following assumptions are now made: 

1. all variances are known; 

2. all variances are assumed equal; 

3. there is no correlation between y^i and y U j, i ^ j and i, j = 1, ... , M and u — l ,..., AT. 
If Assumption 2 alone is violated, appropriate weights could be added in the object. If, in 

addition to Assumption 2, Assumption 3 is also violated, besides adding weights, the cross terms 
(accounting for correlations) should also be included into the objective function. 

Iii the Bayesian framework, the likelihood function (based on observations Y) for 0 and £ can 
be expressed as 


f(<?,£|r)ocp(Ej£,0) 

ocfEj-^exp 



where 

n 

Sij (0i , 0j ) — j 

u=l 

n 

= y [Vui ~ Vi (fui) 0i)] [i luj ~ Vj (£uji Qj)\ 

U=1 

~ * &■> J == 7 TI. 

Box and Tiao [7] used the Bayesian argument setting the prior distribution to be 

p(0,E) =p(%(£) oc [£|-^+^ 2 

based on the assumption of independence between 0 and £. The joint posterior distribution of 
(0, £) is then 


p(0,£|r)ocp(0,£)i(0,£) 


|-(JV+M+l)/2. 


cc |Ej ‘ exp 


It then follows that 

p(0, £ -1 | y) ex |£| < ' JV ~ Af-1 ^ 2 exp 


--trS- 1 ^) 


— trE" 1 ^) 


—oo < 0 < oo, £ > 0, 


—oo < 6 < oo, £ > 0. 

L J 

Since Wishart distribution W m (B~ x ,q) has its probability density function 

W m (B-^q) — p{Z) = fcl^l^-^exp (-\txZB^ , Z > 0, 


where Z is distributed as Wishart with q degrees of freedom and parameter matrix B -1 , it can 
be shown that 



dZ _ j^j|-l/2{9+m-l)2l/2m(g+m-l)p m 
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where 


r ^-[ r G)] n, r ( 6+ —) 


Thus, the posterior density can be integrated over 2 to give 

P(0|S)<xl£ T Er JV/2 . 


b > 


1 — a 
2 ’ 


As a result, the parameter estimates are to be chosen so as to minimize 

d{9) = |£ t E| . (3.2) 

The Bayesian determinant criterion d{8) involves constraints on the number of observations (A r ), 
the number of responses (M), and the number of parameters (P). We need to have N > P and 
N < M in order to carry out the necessary computations; see [8]. 


3.2. Results 


3.2.1. The first-stage analysis 

The first-stage involves a two-compartment model which forms a preliminary work for the 
estimation of parameters of a four-compartment model to be considered later. 

Since the monthly amount of urine cadmium is linearly related to kidney cadmium, the rate 
of change of expected kidney cadmium %(£) at time £ can be represented by 

If kidney cadmium at time £ = 0 is %.(0) and A* 7 ;* is the urinary output rate, the joint model 
for the urine and kidney cadmium is given by 

Vk(t) = ~ (1 - exp (—A*,*)) + 77 / 5 ( 0 ) exp (-A*t), 

=Rk(\- exp (—Afc£)) + \kVk(0) exp (-A /,£), 

where f)k and tj„ denote the expected responses of kidney and urine cadmium, respectively, A*, is 
the rate from kidney cadmium to urine cadmium, and Rk is the monthly intake of cadmium from 
the environment into kidney. The diagram for the two-compartment model of kidney cadmium 
and urine cadmium is presented in Figure 2. 


i 

Kidneys 

i 

Urine 


Figure 2. Diagram for the two-compartment model {Cd_u and CdJc). 

In the least-squares approach based on (3.1), M = 2 and 9 is a 2 x 1 vector of unknown 
parameters, which is [A/,,.i?fcJ T . During the parameter estimation process, measurement of kidney 
cadmium at time 0, which is 12/1983, is added as an unknown parameter. The unit of monthly 
urine cadmium is converted by multiplying by 53 on the assumption that daily output of urine 
is 1-75 g. The results from the nonlinear least-squares estimation show that the range of estimates 
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of At is from 0.0047 to 0.0537, which reveals that it is essential to include the input rate into the 
model. Furthermore, the results indicate that the standard deviations for most subjects are quite 
large, even though the estimates of kidney cadmium at time 0 and prediction of kidney cadmium 
at time 76 are rather close to the measured values. Since the observed values of kidney cadmium 
are quite dominant in the data so that the second assumption of the least-squares approach is 
violated, weights were added into the objective function. 

In order to check the sensitivity of the weights used in the two-compartment model, several 
different ways were used to add weights into the model, based on prior experience with regard to 
the measurement error. 

After weights were added, the residual sum of squares got dramatically reduced. It was in¬ 
teresting to observe that the estimate of A & was not really sensitive to the different methods of 
weights used. In Subject 7, the estimate of A* was the smallest in the model with every weight. 
The largest estimate of Ajt was located either in Subject 11 or 12 among the different weights. 
Table 2 shows one of the results from the weighted least-squares approach. 


Table 2. Results of the nonlinear least-squares for the two-compartment model (Cd_k 
and Cd-u) based on weight, 


Subject 


i 

Cd.k(O) 

Cd k(0) 
(S.) 

Cd_k(76} 

Cd k(76) 
(B.) 

SSQ 

1 

0.007 

40.09 

2303.9 

2500 

5069.61 

4000 

6.933 

2 

0.0167 

262.34 

20965.8 

21000 

23079 

24000 

3.761 

3 

0.016 

254.41 

.. ...... 

13795.7 

12000 

20360.6 

24000 

5.454 

4 

0.0069 

193.22 

15846.9 

14000 

30187.6 

34000 

8.449 

5 

0.0124 

141.13 

3786.7 

2500 

9897.49 

12000 

10.743 

6 

0.014 

345.5 

21041.6 

21000 

30684.4 


5.506 

7 

0.0048 

94.26 

2472.3 

2500 

9435.7 


4.679 

S 

0.0116 

204.19 

15012.9 

15000 

22747.3 


3.230 

9 

0.0220 

425.65 

15253.9 

15000 

21481.3 


2.206 

10 

0.0130 

286.69 

68846.2 

78000 

65108 

24000 

119.27 

11 

0.0271 

1028.3 

31154.9 

31000 

41051.3 

42000 

3.339 

12 

0.0495 

675.8 

10966.0 

10000 

13844.9 


3.417 

13 

0.0150 

681.9 

24928.9 

16000 

46866.5 

66000 

53,046 

14 

0.0236 

630.5 

2395S.0 

24000 

30242.9 


2.258 


S. refers to single. 


In the Bayesian approach, some variables of kidney cadmium were simulated to satisfy the 
Bayesian constraint such that N > P and N > M. It is assumed that the measurement errors of 
kidney fcadmium are normally distributed. The generated variable of kidney cadmium at time t 
is provided by 

Cd-k(t) = nu + n, 

where n was generated to be normally distributed with mean zero and standard deviation one, 
the initial time to is set at 0 and the unit of time is month, fi is the prediction of kidney cadmium 
at time t, derived from the measurements of kidney cadmium at time t = 0 and t = 76, and cr is 
the standard deviation of the measurement error in kidney cadmium, which is assumed to be 
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6000 /ig. Therefore, the simulation is based on the measurements of kidney cadmium at time 0 
and 76, standard normal variables n, the standard deviation of measurement error in kidney 
cadmium cr, and the prediction of kidney cadmium fi. 

Table 3 shows the results from this Bayesian approach, in which the estimated A* for each 
subject falls in a reasonable range that is not too different from the results of the nonlinear 
least-squares. The variation of the estimate of Rk is quite large under both approaches. One 
of the reasons for this may be that the input rate is not directly related to urine cadmium but 
is connected directly to blood cadmium, and yet blood cadmium is not included in the two- 
compartment model. 

Table 3. Results of the two-compartment model from the Bayesian approach with 
the combination of simulated and real data. 


Subject 

Ak 

Rk 

Cdk(b) 

Cd k(0) 

(s.) 

CWk(70) 

Cd k(76) 
(B.) 

O.P.V. 

1 

0.006 

54.50 

2759.63 

2500 

6824.4 

4000 

20.49 

2 

0.018 

276.50 

19169 

21000 

21211 

24000 

38.82 

3 

0.019 

235.82 

11520 

12000 

14920 

24000 

65.22 

4 

0.008 

170.22 

14754 

14000 

25759 

34000 

77.76 

5 

0.015 

105.64 

36780 

2500 

28316 

12000 

280.93 

10 

0.012 

516.62 

71838 

78000 

83474 

24000 

2009.22 

11 

0.027 

990.71 

31958 

31000 

40191 

42000 

21.4776 

13 

0.027 

275.53 

19518 

16000 

13909 

66000 

525.02 


In order to obtain more precise estimates of the rate, the blood and liver cadmium are included 
in the model in our second-stage analysis. From Table 3, it is noted that the estimated Rk with 
low u.p 2 is lower than that with high u~/3 2l which suggests the region of starting values for some 
unknown parameters, at the next step of estimation in the four-compartment model. 

3,2.2. The second-stage analysis 

In the four-compartment model, there are six unknown parameters Cj,[, Cfo, Cbk, Ckb, C* u , 
and R. The required optimization is very complicated in this case of the four-compartment model, 
since the number of responses and the number of unknown parameters are rather large. One of 
the most important things to ensure a successful nonlinear analysis is to obtain good starting 
values for the unknown parameters, before implementing the estimation procedure. There are 
two techniques used for determining the starting values of unknown parameters in this four- 
compartment model, viz. a graphical analysis and a grid search. 

With the graphical analysis, the meaning of each parameter can be interpreted graphically 
to describe the behavior of the expectation functions. There are three references, the half-life 
times of liver and kidney cadmium and blood compartment with a rapid turnover from previous 
results, the portions of body burden of kidney and liver cadmium from previous results, and 
the estimates of excretion rate and cadmium intake from the two-compartment model discussed 
earlier, forming the basis of interpreting the unknown parameters [Cbi, Cm, Cbk, Ckb, Cku, R}'- 

The steps involved in the graphical analysis are as follows. 

1. Solve the expectation functions numerically after fixing the parameters within reasonable 
ranges based on the meaning of the parameters just mentioned above; Matlab software 
package was used for obtaining the numerical solutions, using Runge-Kutta-Felilberg it¬ 
eration method. 
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2. Draw the sketches of both fitted values derived from the numerical solution and the ob¬ 
served values of blood cadmium against time, liver cadmium against time, kidney cadmium 
against time, and urine cadmium against time, simultaneously, and visualize them. 

3. Change the values of fixed parameters and repeat steps (l)-(3) until the fitted and ob¬ 
served values of blood cadmium, liver cadmium, kidney cadmium, and urine cadmium are 
reasonably close. 

With the grid search, a combination of all resources for the graphical analysis and the results of 
the graphical analysis were used as references for the range of starting values of the parameters. 
In attempting to find proper starting values for the grid search, there are many steps taken. In 
the beginning, the increments are quite large for a rough search; then the increments become 
smaller and smaller for a finer search each time. The program written in C is intended not to 
give final results automatically, but to provide middle results step by step. During the search, 
it can be visualized on screen in order to give an idea of which part of the range is most likely 
to be selected. Moreover, the program can be paused at any time in the middle of searching 
if it is necessary, for deciding about adjusting the increments or even changing the range of 
parameters. The results of the grid search for the starting values of the parameters are shown 
in Table 4. Although this type of searching is very time consuming, it provides a general idea 
of where the location of the local minimum might be obtained within a reasonable range of the 
parameters. Furthermore, it most certainly accelerates the convergence of the solution of the 
nonlinear, optimization problem. 

Table 4. Results of the starting values for the parameters obtained from grid search 
for the four-comportment model. 




6 lb 

&bk 

c kb 

Cku 

ft 

O.F.V. 

Subject 1 

22 

.002 

7 

.017 

.007 

414 

30.5 

Subject 2 

7.8 

.0033 

13.1 

.0045 

.017 

330 

17.83 

Subject 3 

17 

.043 

18.4 

.0092 

.018 

141 

25 

Subject 4 

13 

.0001 

12 

.0038 

.008 

512.5 

19.96 

Subject S 

0.2 

.0019 

14.2 

.0001 

.0105 

121.1 

36.6327 

Subject 10 

17 

.007 

1.8 

.0004 

.0116 

530 

154.81 

Subject 11 

9.6 

.0059 

19.2 

.015 

.027 

1265 

9.3449 

Subject 13 

8.4 

.0013 

11.2 

.0035 

.03 

840.7 

100.9133 


Table 6 shows the results of the four-compartment model based on the nonlinear least-squares 
estimation, while Table 6 shows the results of the four-compartment model from the Bayesian 
approach based on simulated and observed data. The steps of generating the simulated data 
are similar to those based on the two-compartment model. Estimates of Cm, Co,, Cj,*, and C^b 
from both approaches provide some information about the blood distribution of cadmium to 
liver and kidney and the cadmium elimination rate from blood either to liver or to kidney, which 
the two-compartment model was unable to provide. By comparing the estimates of Cku hi the 
four-compartment model and of A* in the two-compartment model, it is noticed that the results 
from the two-compartment model and the four-compartment model remain close. However, the 
estimate of Rk in the two-compartment model differs greatly from the estimate of R in the four- 
compartment model. The convergence of the optimization procedure goes smoothly due to the 
fact that the grid search provides good starting values for the unknown parameters. 

All subjects except 13 show that the fitted data from the four-compartment model agree 
reasonably well with the observed data. u.p 2 for Subject 13 is very high, which indicates that 
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Table 5. Results of the nonlinear least-squares estimation method for the four- 
compartment model. 


Subject 

C bl 

Cib 

Cbk 

Okb 

Ckn 

ft 

O.F.V. 

1 

21.98 

.00001 

6.993 

.0164 

,0067 

414.0 

27.7 

2 

8.68 

.0039 

21,32 

.013 

.0169 

325.03 

16.7581 

3 

14.30 

.0568 

15.70 

.0037 

.0177 

141.00 

24.20 

4 

13.39 

.00001 

12.24 

.0038 

.0079 

510.09 

19.86 

5 

6.74 

.0705 

10.82 

.0001 

.0107 

62.10 

12.0684 

10 

17.25 

.0075 

1.15 

.0002 

.0116 

530.0 

154.11 

11 

12.02 

.008 

17.98 

.013 

.0273 

1259.9 

9.22 

13 

9.15 

.0021 

11.55 

.0045 

.0304 

837.21 

100.82 


Table 6. Results for the four-compartment model from the Bayesian approach based 
on simulated and observed data. 


Subject 

Cu 

<%b 

Cbk 

Ckb 

Cku 

R 

O.F.V. 

2 

10.33 

0.0054 

19.67 

0.0116 

0.0171 

294.82 

3216.74 

4 

14.05 

0.000004 

11.89 

0.0033 

0.0078 

523,5 

5668.0 

11 

11.46 

0,0085 

18.54 

0.0166 

0.0253 

1174.09 

79.44 


the kidney functions for this subject may be damaged and some irregular and unpredicted pattern 
could have taken place in this case. Figures 3~5 are examples showing good agreement of the 
observed data with the fitted four-compartment model for Subjects 2, 4, and 11. 

3.3. Discussions and Suggestions 

The large estimates of Cfo and Cbk confirm that there is no cadmium accumulation in blood 
with rapid turnover of cadmium concluded from previous results. 

A combination of Ckb and Cku reflects the half-life time in the kidney. Absence of feces cadmium 
variable might cause the estimate of the half-life time in the kidney to be small. It is probable 
that missing observations of feces cadmium leads to the overestimation of Ckn¬ 
it is noticed from the results that estimates of Cib, Cbi, Cbk> Cjtu, and R might be affected by 
the absence of observations of tissue cadmium and feces cadmium. The range of estimates of Ca 
indicates the half-life time in liver. The estimate of Ca, seems small, which might be due to no 
tissue cadmium variable being involved in the model. It might affect the estimate of Cib since 
cadmium may be transported back from tissue to liver via blood, and missing observations of tis¬ 
sue cadmium might lead to both Cu and Cbk being overestimated and Qj, being underestimated. 
Unfortunately, there is no entrance to measure the tissue cadmium yet. Therefore, when the 
half-life time of blood, liver, and kidney cadmium are all derived from each of these estimates, 
the influence of the absent tissue cadmium should be taken into account. 

Since the tissue cadmium is about 31% of the total body burden, it should be included in the 
model when the number of measurements of the blood cadmium and others increase to a certain 
extent, in spite of unavailable observation on the tissue cadmium. The half-life time of tissue 
cadmium and the relationship with the other compartments can be estimated through blood 
cadmium and other observations without tissue cadmium, based on the Bayesian approach, in a 
similar way as shown in Table 7. Table 7 gives the results of the four-compartment model from 
the Bayesian approach derived just from blood cadmium observations. It should be noticed that 
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Figure 3. Plot of the observed data for Subject 2 with the fitted four-compartment 
model. 



50000 
400005' 
4- 30000*: 
20000J 
10000 
0 


Figure 4. Plot of the observed data for Subject 4 with the fitted four-compartment 
model. 


Table 7. Results For the four-compartment model from the Bayesian approach based 
only on blood cadmium observations. 


Subject 

On 

On 

Oik 

Okt 

Oku 

R 

O.F.V. 

1 

22.0 

.000003 

7.00 

.0177 

.0071 

414.0 

4.24 

2 

7.8 

.0039 

13.10 

.0070 

.0531 

330.0 

6.29 

3 

14.3 

.062 

15.70 

.0064 

.0375 

140.90 

7.48 

4 

14.0 

.000006 

12.00 

.0054 

.0105 

500 

2.04 

5 

6.70 

.0704 

11.00 

.0002 

.000008 

60 

1.19 

10 

16.98 

.006 

1.79 

0 

.007 

739 

6.88 

11 

9.40 

.013 

16.79 

.00001 

.816 

1265.0 

1.77 

13 

8.40 

.000004 

11,2 

.0051 

.018 

840.70 

9,93 


the results derived from blood cadmium alone are based on good starting values provided by the 
graphical analysis and the grid search based on all responses. 

The estimate of R is quite reasonable. It is found from the results of the four-compartment 
model that subjects with high /? 2 -microglobulin have higher estimates of intake cadmium (R) 
than those with low /3 z -microglobuiin. The interpretation for this could be that the environment 
is not only the source for cadmium intake, but that the cadmium released from other tissues 
might also contribute to the intake of cadmium. 

The differences of other estimates like <?(,[, Cu,, Cm, Ckt. and C^u between the two subgroups 
are not obvious. It is not conclusive because of a small and insufficient dataset and the uncertainty 
in observations; it perhaps implies that the distribution of cadmium in the body depending on 
the absorption and the metabolism of cadmium in the body varies from individual to individual, 
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Figure 5. Plot of the observed data for Subject 11 with the fitted four-compartment 

model. 

all that it is not statistically different between a worker with abnormal kidney condition and one 
with normal kidney condition. 

In order to formulate further refined models and develop more powerful statistical analysis for 
the model* as the cadmium research still continues on actively, there are four suggestions for 
cadmium researchers. First, the number of subjects and the number of measurements should be 
increased in order to carry out appropriate statistical analysis effectively. Second, since it is more 
economical to collect blood and urine cadmium samples from workers exposed to cadmium than 
liver and kidney cadmium in vivo measurements, blood and urine samples should be gathered 
more often for preliminary analysis. A series of blood cadmium variables can be used for rough 
estimation of Cbi, Cat, and it by a four-compartment model. Even though the other three 
parameters, Cib, Ckbi and Cfc u , can also be estimated through blood cadmium, the estimation 
might be poor for Cjt, and C*{. (with very small values) and for Ck u (with no direct connection 
to blood cadmium), but they can Still be a reference for researchers. In the meantime, a series of 
urine cadmium variables and a few kidney cadmium variables can be used to estimate R and C* u 
by a two-compartment model. In addition, tissue cadmium variable can be incorporated into the 
model in order to estimate some rates relating to the tissue cadmium. Third, it may be possible 
in the future to investigate the cadmium metabolism in the human body in stratified subgroups 
with kidney condition when sufficient data are available. Moreover, nonlinear segmented models 
according to different exposure levels in different time periods can be formulated. Consequently, 
several statistical analyses may be needed based on different time intervals. The model should 
be formulated on the data with shorter time intervals at first, and then designed for the data 
with longer and longer time intervals. Finally, after the completion of in vivo measurements, 
entire analysis can be done using a four-compartment model through the Bayesian approach. 
It is also possible to analyze the cadmium metabolism by a five-compartment model including 
the tissue cadmium, provided sufficient in vivo measurements of liver and kidney cadmium are 
available. Since different responses in the same case would have been collected from the same 
body, correlation will exist between them. Therefore, we strongly recommend the use of the 
Bayesian determinant minimization method for the cadmium analysis whenever sets of such 
multiresponse data are available. 
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